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ABSTRACT 


Presented  in  this  report  are  finite-difference  (forward)  heat 
transfer  equations  applicable  to  transient,  radial  heat  flow  in  spheres 
and  cylinders  and  applicable  to  transient,  one -dimensional  heat  flow  in 
flat  plates.  By  very  minor  manipulations,  a  single  heat  transfer  procedure 
can  easily  be  utilized  to  determine  transient  heat  flow  in  all  three  basic 
structure  configurations. 

For  each  skin  configuration  the.  accuracy  of  the  finite-difference 
procedure,  compared  with  exact  analytical  methods,  depends  on  optimum 
selection  of  the  calculation  time  increment  and  the  incremental  distance 
between  temperature  nodes  in  relation  to  the  material  thermal  properties 
and  on  the  closeness  of  the  approximate  temperature  gradients  to  the  true 
gradients . 
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Stction  I.  INTRODUCTION 


The  U.  S.  Army  has  had  great  success  using  finite- difference 
(forward)  techniques  to  accurately  determine  transient  heat  transfer 
in  homogeneous  and  composite  wall  missile  structures.  Heating  effects 
on  most  Army  missile  structures  have  been  analyzed  by  use  of  finite- 
difference  methods  applicable  to  transient,  one- dimensional  heat  flow 
in  flat  plates.  This  flat-plate  technique  is  generally  quite  adequate 
for  structure  analyses  of  missiles  such  as  JUPITER,  PERSHING,  REDSTONE, 
and  SERGEANT  because  most  of  the  airframe  exposed  to  inflight  aero¬ 
dynamic  heating  has  a  relatively  large  radius  of  curvature  in  relation 
to  the  skin  thickness.  Thus  the  curved  structure  actually  approached 
a  flat  plate  in  relation  to  heat  storage  and  heat  flow. 

Contrarily,  a  flat-plate  solution  may  not  suffice  for  analysis  of 
small  hemispherical  tips,  small  hemicylindrical  leading  edges,  small 
caliber  vehicles,  large  or  small  cylindrical  and  spherical  solid 
propellant  motors  exposed  to  environmental  thermal  shocks,  and  many 
other  missile  components  and  sections.  Thus,  accurate,  flexible 
procedures  for  at  least  three  basic  structure  configurations  are 
necessary  for  proper  thermal  analyzation  of  present  and  future  Army 
missile  systems. 

One  purpose  of  this  report  is  to  focus  attention  on  the  similarities 
existing  in  finite-difference  heat  transfer  equations  for  radial  heat 
flow  in  spheres  and  cylinders  and  one-dimensional  heat  flow  in  flat 
plates.  Another  purpose  of  this  report  is  to  present  one  general  heat 
transfer  procedure  which  can  be  used  to  determine  transient  heat  transfer 
effects  in  three  basic  configurations,  i.e.,  sphere,  cylinder,  and  flat 
plate . 
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Section  II.  FLAT-PLATE,  CYLINDER,  AND  SPHERE  PROCEDURES 


1.  Background 

Transient,  finite-difference  (forward)  heat  transfer  methods 
have  occupied  a  vital  role  in  the  development  of  several  of  the  Army's 
most  prominent,  reliable  missile  systems.  During  the  historical  heat 
protection  material  development  program  for  the  JUPITER  nose  cone, 
flat-plate,  finite-difference,  heat  transfer  methods  were  used  success¬ 
fully  in  analyses  of  reentry  simulation  test  results1  and  in  analyses 
of  inflight  aerodynamic  heating  effects.2'3  This  same  heat  transfer 
method,  in  conjunction  with  a  simultaneous  ablation  routine,4  was  used 
exclusively  to  determine  the  thermal  protection  requirements  for  the 
PERSHING  reentry  body.  Later,  calculated  inflight  temperatures  were 
found  to  be  in  good  agreement  with  temperatures  measured  on  the  reentry 
body  of  PERSHING  missile  308. 5  Agreement  of  finite-difference  results 
with  exact  solutions  is  discussed  briefly  in  Paragraph  6  of  this 
section. 

With  the  advent  of  large  solid-propellant  power  plants  for  use  in 
severe  environments,  as  well  as  increased  demand  for  hypersonic,  low- 
altitude  flights  by  small  vehicles  or  large  vehicles  having  small  lead¬ 
ing  edges  or  tips,  the  Army  has  become  increasingly  involved  with  heat 
flow  problems  which  are  not  adequately  solvable  by  flat-plate  techniques. 
Recently,  therefore,  transient,  radial  heat  transfer  procedures  for 
cylinders  and  spheres,  in  finite-difference  form,  have  been  derived. 

Using  tu_  TSula’i.  tut  5  CyM'i.ict ,  gvcd  agTetiii  Ifetwtr?!. 

measured  and  calculated  temperatures  of  a  PERSHING  XB-1  motor  subjected 
to  severe  thermal  shock  in  a  cooling  chamber  were  obtained.6 

Procedures  for  determining  radial  heat  flow  in  cylinders  and 
spheres  have  many  applications  in  missile  design  and  analysis,  since  in 
numerous  cases  the  assumption  of  only  radial  flow  is  sufficiently 
accurate.  For  instance,  transient,  radial  heat  flow  in  a  small,  solid 
hemispherical  tip  may  be  adequate  for  materials  with  low  thermal 
diffusivities  (reinforced  plastics,  for  example)  even  when  the  heat 
input  varies  with  external  surface  station.  Transient,  radial  heat 
flow  at  the  midsection  of  a  cylindrical  solid-propellant  motor  having 
a  large  L/D  ratio  may  be  adequate  in  most  cases  because  the  present 
solid  propellants  have  small  thermal  diffusivities  and  end  effects 
do  not  influence  the  motor  midsection  for  many  hours.6 

The  conduction  equations  presented  for  each  of  the  three  basic 
structure  configurations  are  applicable  to  any  number  of  material 
layers;  however,  only  two  layers  are  illustrated  because  this  is  the 
maximum  number  of  materials  necessary  for  the  application  of  all 
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equations  for  thermally  thick  skins.  Equations  for  thermally  thin 
skins  and  thick- thin  combinations  have  been  derived  for  flat  plates, 
cylinders,  and  spheres  to  account  for  all  practical  material  arrange¬ 
ments.  These  equations  are  not  presented  due  to  similarity  of  the 
derivations  to  those  of  the  thick- thick  wall  cases  discussed. 

Many  statements  made  throughout  the  discussion  of  the  flat-plate 
procedure  are  also  applicable  to  the  other  two  procedures,  and  are 
not  repeated  to  preclude  repetition.  Since  the  average  area  of  a  sphere 
across  which  heat  flows  is  not  a  linear  function  of  the  radius,  the 
derivation  of  average  areas  for  the  spherical  configuration  is  shown  in 
Appendix  A.  Shape  factors  have  been  purposely  omitted  from  the  radia¬ 
tion  terms  throughout  the  report. 


2.  Transient,  OnrDimtniional  Heat  Transfer  for  Flat  Plates 

One-dimensional,  flat-plate  heat  transfer  in  a  homogeneous 
material  can  be  determined  by  solving  heat  balance  equations  at  the 
exposed  surface,  unexposed  surface,  interior  nodes,  and  interfaces. 

The  forward  finite  difference  method  was  used.  It  was  assumed  that 
the  incremental  thickness  (r)  can  be  selected  sufficiently  small  to 
give  accurate  temperature  gradients  between  adjacent  nodes  and  that  the 
incremental  time  (At)  is  small  enough  to  neglect  any  effect  on  regions 
more  than  one  r  from  the  node  in  question. 
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Figure  1 


a .  Exposed  Surface 
From  Figure  1  the 

^convQ  +  ^solar  "  ^radQ 


heat  balance  at  the  exposed  surface  is 
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where 


qconv0"  hoA  ^Teff0  “  Ti) 
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1 


'a 


!§ A(T* : Ti> 
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The  area,  A,  is  uniform  for  one- dimens ional ,  flat-plate  heat  transfer. 
Rewriting  Equation  (1)  we  have 

kaO*  *  Tx)  _  Ta  (T {  -  TX) 
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Multiply  by  ,  let  0  =  ,  let  Y  -  ~  ,  and  solve  for  Tx 

T{  -  Ti(l  -  2Yah0  -  2@a)  +  20aTa  +  2Yah0Teff^  +  2Ya«0a(T*o  -  T?)  + 
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The  coefficient  of  Tx  in  Equation  (3)  must  be  >  0  otherwise 
depend  on  Tx  in  the  negative  sense  wMflt  i*  tot 
set  the  coefficient  of  Tx  equal  to  zero  and  solve  for  S  . 
Equation  (2)  that  Y  ■  0<T/k).  a 
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Tlifcteiort, 
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EquatioiWR^f^  ^  Crltlcal  (0a>  and  the  upper  Unit  of  the 
l  in  Equation  (8)  for  convergence.  It  is  readily  seen  that  if  0  in 

Equation  (8)  is  a  maximum  of  0.5  (to  make  the  coefficient  of  T  -  0) 

C°®fficjent  of  Ti  in  Equation  O)  is  negative  because^Yh  is 
positive.  Therefore,  0  must  be  equal  to  or  less  than  \ 
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(4) 


The  maximum  h0  expected  or  calculated  for  the  specific  case  for  analysis 
should  be  used  in  determining  the  critical  0. 


b.  Interface  Equation 


the  energy  balance  is 


Q  -  fl  an 

cond  ^cond  ^stored 


At  the  interface  between  material  "A"  and  "B"  (T  Figure  1) 

1  an/^a  4  a  ®  ' 


n-i-^n  n-*n+i 
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Rearrange  and  solve  for  T' 
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2k  At 
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c .  Internal  Nodes 


The  heat  balance  at  any  interior  point  of  a  homogeneous 
wall  (Figure  1)  may  be  written 
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Solving  for  T^+l  with  9  -  gives 
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d.  Backside  Surface 


written  as 


The  energy  balance  at  the  backside  surface,  T.  ,  may  be 


‘stored 


bs-l  —  bs 


kb^Ib»-l  *  Tba* 
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Rearrange  and  solve  for  T^g 
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To  find  the  limiting  0  for  the  backside  material,  equate  the  coefficient 
°f  Tbg  (Equation  10)  to  zero  and  solve  for  0b> 


3.  Transient,  Radial  Heat  Transfer  for  Cylinders 


a.  External  Surface  Equation 


Consider  a  cylindrical  segment  heated  as  shown  in  Figure  2. 
From  the  energy  balance  at  the  peripheral  surface 
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rearrange,  and  solve  for  . 
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Equating  coefficient  of  T:  to  zero,  since  cannot  depend  on  the 
previous  Tj  in  the  negative  sense; 
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Since  Y  =  3  7 
k 


the  critical  or  limiting  g  at  the  external  surface  is 


(14) 


For  R  — 00  and  T  finite  or  x/R«  1,  Equation  (14)  approaches  the  one- 
dimensional  flat  plate  case (Equation  4) . 


b.  Interface  Equation  (Figure  2) 

At  point  Tn  (interface  between  Material  "A"  and  "B")  the 
energy  balance  becomes 
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(15) 
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Figure  2 
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solving  for 


T' 

n 


2ka4t(R  -lTn  +^) 

T"  +  Ta(R  -  ET„)(paCaTa  +  pbCbTb)  <Tn-l 

2kbZjt(R  ■  £■,„  - 

Tb(R  -ITn)^PacaTa  +  PbcbTb>  (Tn+1  '  Tn> 


(16) 


c.  Interior  Equation  (Figure  2) 

At  point  Tn^  (Material  "B")  the  energy  balance  becomes 
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Solving  for  the  critical  £  at  the  inner  surface  yields 


4  J  L 


For  small  Tb's  and  small  ZTbg's,  the  one  dimensional  flat  plate  case 
(Equation  11)  ic  approached  as  R  becomes  relatively  large. 
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e.  Observations 


For  finite  r's  and  R-^oo,  cylindrical  temperature  Equations 
(13),  (16),  (18),  and  (20)  revert  to  flat-plate  temperature  Equations 
(3),  (6) ,  (8),  and  (10),  respectively. 


4.  Transient,  Radial  Haat  Tran*l#r  lor  Sphoroi 


Figure  3 


a .  External  Surface  Equation 

Consider  a  spherical  segment  heated  as  shown  above. 
From  the  energy  balance  at  the  pheripheral  surface 

qconv0  +  qsolar  "  qrad0  "  q^ond  “  ^stored 

l-*2  1 
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Note:  See  Appendix  A  for  derivation  of  average  areas. 
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Set  the  coefficient  of  Ti  equal  to  zero  and  solve  for  the  limiting 
3  when  h  is  a  maximum. 
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For  finite  Ta's  Equation  (24)  approaches  the  flat-plate  Equation  (4)  as 
R  becomes  large. 

b.  Interface  Equation 


balance  is 


At  point  Tn  (interface  of  material  "A"  and  "B")  the  energy 
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Solve  for  T' 
n 


At  point  Tn+1  (Material  "B")  the  energy  balance  is 

^cond  ^cond  "  ^stored 
n-*  n-K  n+i""*’ n+8  n+i 


or 
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d.  Backside  Surface  Equation 


At  the  backside  surface  (T.  )  the  heat  balance  is 
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Equate  the  coefficient  of  Tbg  to  zero  and  solve  for  the  limiting  0b 
when  h^  is  a  maximum. 


As  with  the  cylindrical  Equation  (21)  this  f3b  approaches  the  flat- 
plate  limits  (Equation  11)  when  rb  is  finite  and  R  is  relatively  large. 

c .  Observations 

For  finite  t's  and  R —  *,  spherical  temperature  Equa¬ 
tions  (23),  (26),  (28),  and  (30)  revert  to  the  flat-plate  temperature 
Equations  (3),  (6),  (8),  and  (10),  respectively. 


5.  General  Equations  Applicable  to  Flat  Plates,  Cylinders,  and  Spheres 

Due  to  striking  similarities  between  finite-.iif ference  heat 
transfer  equations  for  flat  plates,  cylinders,  and  spheres,  general 
equations  adaptable  to  all  three  configurations  are  derived  and 
presented  as  Equations  (32)  through  (37).  The  coefficients  G  and  m  and 
the  exponent  m  are  inserted  to  perform  manipulations  resulting  in 
equations  applicable  to  the  desired  skin  configuration.  The  following 
table  gives  values  assigned  to  G  and  m. 


Skin  Configuration 

G 

m 

Flat  Plate 

1 

1 

Cylinder 

2 

1 

Sohere 

2 

2 
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6.  Finite-Difference  Results  and  Procedure  Selection 


Calculated  results  from  the  three  finite-difference  heat 
transfer  procedures  presented  in  this  report  have  been  compared  with 
the  exact  solutions.7*8  Excellent  agreement,  within  0.3  percent, 
was  obtained  in  each  case  for  a  material  initially  at  uniform  tempera¬ 
ture  (T0)  and  with  one  of  its  surfaces  maintained  at  a  constant  ' 
temperature  different  from  T0  after  time  zero.  The  accuracy ‘of  any 
finite-difference  heat  flow  equation  depends  on  the  proper  selection  of 
the  time  interval  (At)  and  distance  increment  (r)  as  well  as  the  close¬ 


ness  of  the  approximate  temperature  gradients, 
gradients . 


Tx 


*  Tx+i 


to  the  true 


For  the  majority  of  missile  structure  analysis,  the  simple  flat- 
plate  procedure  may  be  sufficiently  accurate.  For  instance,  examination 
of  the  cylindrical  and  spherical  conduction  equations  (Equations  20  and 
30)  at  the  inside  surface  (where  the  maximum  deviation  from  a  flat 
plate  obviously  will  occur)  shows  that  a  ratio  of  R/lrbs  (radius /wall 
thickness)  of  100  gives  results  approximating,  within  0.6  percent, 
the  flat-plate  solution  (Equation  10).  Thus,  for  adequate  engineering 
analyses  and  structural  design,  the  simple  one-dimensional  flat-plate 
procedure  is  advantageous  in  many  investigations  because  of  the 
required  computer  time. 
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Section  III.  CONCLUSIONS 


One  general  finite-difference  heat-transfer  procedure  can  be 
used  with  minor  manipulations,  to  calculate  transient,  radial  heat  flow 
in  spheres  and  cylinders  and  to  calculate  one- dimensional  heat  flow  in 
flat  plates. 

Excellent  agreement  between  finite- difference  calculations  and 
exact  analytical  results  is  obtained  when  proper  time  and  distance 
increments  are  chosen  in  relation  to  material  thermal  properties. 

For  hollow  or  solid  spheres  and  cylinders  composed  of  one 
homogeneous  materialithe  limiting  modulus, 0  ,  may  be  determined  by 
parameters  at  the  internal  surface  even  though  the  external  surface  ma, 
be  the  only  surface  exposed  to  convective  heating. 

Comparisons  of  finite-difference  calculations  with  exact  analytical 
results  for  several  special  heat  transfer  cases  are  necessary  for  final 
determination  of  the  accuracy  of  the  finite-difference  methods.  Results 
of  these  additional  comparisons  will  be  published  in  a  future  document. 
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Appendix  A 

DERIVATION  OF  AVERAGE  AREA  EQUATIONS  FOR  SPHERICAL  CONDUCTION 

Since  the  surface  area  subtended  by  a  solid  angle  on  a  sphere  is 
not  a  linear  function  of  the  radius  (R)  the  average  area  (7)  through 
which  heat  flows  and  the  average  area  for  determining  heat  storage 
volume  (  T J  )  must  be  found. 

1.  External  Surface 

The  average  area  from  1  to  2  (Figure  3)  is  found  from 

rR 

cp  R2dR 

*R-T  17  T  \ 

*-—r-  -  9  R-f)  +fj  <38> 

a  L'  '  J 

To  find  the  average  area  for  determining  the  heat  storage  volume  from  R  to 
T 

R  -  the  following  equation  is  used. 


Results  of  the  integration  in  Equations  (38)  and  (39)  are  reflected 
in  Equation  (22)  of  the  text. 

2.  Interface 


The  average  area  from  n-l  to  n  (Figure  3)  is 


and  the  average  area  from  n  to  n+i  is 


J,R-  Er„ 
R3  dR 

R-  Et  -T. 


♦('A'f)  +"I2 


The  average  area  for  finding  the  heat  storage  volume  at  an  interface  (n) 

is  _ 

/»R-  Et+t  /2 


)  I  RSdR 

Jr^Lt  -t  /2 
n  b 


T  +  T, 
a  b 


(R  *  ^Tn  +  ~t)  +  ~48 


.2  _  3 

*-^-r)  +-48 


Results  of  the  integration  of  Equations  (40),  (41),  and  (42)  are  used  in 
Equations  (26). 

3.  Interior 

The  average  area  through  which  heat  flows  from  n+1  to  n+2 
is  given  by  Equation  (43).  The  A  between  n  and  n+ i  is  found  by  using 
+  r  / 2  in  Equation  (43). 


JiR'  “Tn+l 
R2dR 

R-  Et  -T 

n+1  b 


-  9lR-ET, 


T  \2  T  2 

-i)  +i- 

n+1  2/  12 


To  obtain  the  heat  storage  volume  at  node  n+ 1  the  following  average 
area  is  used. 


R-  EWTh/2 


/•“■  “  n+l  b' 
9  I  RadR 

*'R  ETn«’Tb/: 


<{MOa4] 


Results  of  the  integration  in  Equations  (43)  and  (44)  are 
inserted  in  the  interior  Equation  (27). 
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4.  Bac ki id*  Surface 

The  average  area  through  which  heat  is  conducted  from  nodes 
bs-i  to  bs  (Figure  3)  is 


<P 


+TV 


(45) 


and  the  area  determining  the  heat  storage  volume  is 


The  results  of  the  integrations  of  Equations  (45)  and  (46)  are  reflected 
in  Equation  (29). 
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Appendix  B 


LIMITS  OF  THE  DIMENSIONLESS  MODULUS,  3  ,  FOR 
FLAT  PLATES,  CYLINDERS,  AND  SPHERES 

In  all  conduction  equations  presented  in  this  report,  the  modulus, 
3,  has  limits  which  are  critical  to  the  stability  of  the  equations. 

This  is  due  to  the  necessity  of  coefficients  of  a  particular  temperature 
point  being  equal  to  or  greater  than  zero  because  it  is  unreasonable 
for  a  temperature  to  depend  on  a  previous  temperature  in  the  negative 
sense.  The  following  equations  show  the  limits  of  3  under  conditions 
where  the  heat  transfer  coefficient  (h)  equals  zero. 

1.  Heat  Transfer  Coefficient  (h)  =  0 
a  .  Flat  Plate 


The  equations  of  3  applicable  to  the  external  and  internal 
surface  (Equations  4  and  11)  reduce  to  the  following  upper  limit 
when  hG  =  0  and  h^  =  0. 

t  i  |  (47) 

b .  Cylinder 

(1)  External  Surface 

The  limiting  3  equation  (Equation  14)  at  the  external 
surfaces  reduces  to  the  following  when  h0  =  0. 


(48) 


which  has  values  approaching  0.5  and  0.75  for  R»  xa  and  Ta  -*■  R 
respectively. 


(2)  Internal  Surface 

The  limiting  3  at  the  internal  surface  or  center  of  a 
cylinder  (Equation  21)  reduces  to  the  following  when  hj[  =  0. 


2 


/R  -  iTba 


+  V2  \ 
+  V4  / 


(49) 


The  values  of  Equation  (49)  approach  0.25  and  0.5  as2n>s  ^ ^  an® 
when  K'»  ^Tbs  f®aPectlvely.  Thus  when  h^-  0  the  limiting  modulus  0  of 
a  cylinder  composed  of  one  homogeneous  material  may  be  determined  by 
Equation  (49)  rather  than  Equation  (48).  The  limiting  0  varies  radically 
from  0.5  only  when  R  -  is  less  than  one  r. 

c.  Sphere 

(1)  External  Surface 

The  critical  0  in  Equation  (24)  reduces  to  the  follow¬ 
ing  when  hQ  =  0. 


The  values  of  0  in  this  equation  approach  0.5  and  0.875  when 
R»ra  and  when  ra— ►R  respectively. 


(2)  Internal  Surface 

The  limiting  0  in  Equation  (31)  reduces  to  the  following 

when  h.  =0. 
l 


(51) 


The  0  in  this  equation  has  values  approaching  0.125  and  0.5  when 
2rbg— ►R  and  when  R»  2rbg  respectively.  Thus  the  overall 

critical  0  for  solid  or  near  solid  spheres  composed  of  one  homogeneous 
material  may  be  determined  by  Equation  (51)  rather  than  Equation  (50). 
The  lower  0  limit  of  0.125  is  approached  rapidly  after  R  -  2rbg  becomes 
less  than  one  r. 


2.  Htot  Transfer  Coefficient  4  9 

To  determine  the  limits  of  0  for  each  of  the  three  configurations 
when  h  f  0  the 0  must  be  determined  for  the  maximum  heat  transfer 
coefficient..  Both  external  and  internal  surfaces  0's  should  be  calculated 
to  determine  upper,  limits  which  are  valuable  in  insuring  stability  or 
convergence^  of  the  equations. 
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and  applicable  to  transient,  one-dimensional  heat  flow  in  flat  plates. 

By  very  minor  manipulations,  a  single  heat  transfer  procedure  can  easily 
be  utilized  to  determine  transient  heat  flow  in  all  three  basic  structure 
configurations. 

For  each  skin  configuration  the  accuracy' of  the  finite-difference 
procedure,  compared  with  exact  analytical  methods,  depends  on  optimum 
selection  of  the  calculation  time  increment  and  the  incremental  distance 
between  temperature  nodes  in  relation  to  the  material  thermal  properties 
and  on  the  closeness  of  the  approximate  temperature  gradients  to  the  true 

gradients. 
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